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We employ the ATLAS search results for events containing jets and large missing 
transverse momentum, corresponding to an integrated luminosity of 1 fb~^, to inves- 
tigate the constrained minimal super symmetric model (CMSSM) with h — r Yukawa 
coupling unification. In this scenario, one of the stops is the next-to-lightest super- 
symmetric particle (NLSP), which co-annihilates with the lightest (LSP) neutralino 
to yield the desired dark matter relic abundance. The NLSP stop, here taken to be 
lighter than the top quark, is slightly 20% — 30%) heavier than the LSP neu- 
tralino, and it primarily decays into the LSP and a charm quark. We find that the 
multi-jets and monojet ATLAS searches are sensitive to this scenario if the stop pair 
production is accompanied by a hard QCD jet. The excluded limit for the NLSP 
stop mass from the ATLAS data can reach 160 GeV in the coannihilation region, 
with mass below 140 GeV essentially excluded. A significant region of the parameter 
space corresponding to large mo values, 8 TeV ^ mo ^ 16 TeV, is excluded by our 
analysis. For LSP neutralino mass ~ 100 GeV, the LHC constraints in some cases 
on the spin-dependent (spin-independent) neutralino-nucleon cross section are sig- 
nificantly more stringent than the current and expected bounds from Xenon, CDMS 
and IceCube. 
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I. Introduction 



Low scale supersymmetry, augmented by an unbroken R-parity, largely overcomes the 
gauge hierarchy problem encountered in the Standard Model (SM) and also provides a com- 
pelling cold dark matter candidate. In the mSUGRA/constrained minimal supersymmetric 
model (CMSSM) as well as in many other realistic models, the lightest neutralino (LSP) 
is stable j^j with a relic density that is compatible with the WMAP dark matter measure- 
ments [sl. However, the small annihilation cross section of a pure bino LSP with mass of 
around 100 GeV does not permit one to easily reproduce the required relic dark matter 
abundance j^j- An interesting scenario which overcomes this conundrum is bino-NLSP stop 
coannihilation. In this case the bino and the NLSP stop are quasi-degenerate in mass, such 
that the ensuing coannihilation processes in the early universe allow one to reproduce the de- 
sired dark matter relic density. Other related scenarios include bino-gluino coannihilation jsj 
and bino-sbottom coannihilation jof. 

Besides the coannihilation scenario, there are other motivations for considering an NLSP 
stop. For instance, successful electroweak baryogenesis seems possible in the presence of a 
light stop and favors M^^,Mh ^ 130 GeV Also, the stop-top loop diagrams contribute 
to the parameter m'jj^, and naturalness in electroweak symmetry breaking in the MSSM 
consequently prefers an upper bound on the lighter stop mass of around a few hundred 
GeV |8|. The larger top Yukawa coupling typically yields third generation squark masses 
that are lighter than the first two generations jol. 

Our study here is inspired by a recent study 10| which showed that SU(5) or SO(IO) 
inspired b — t Yukawa unification is compatible within the CMSSM framework with the 
WMAP dark matter bounds only if there exists NLSP stop - LSP neutralino coannihilation. 
The universal scalar parameter mo in this scenario turns out to be of order 5 — 20 TeV, with 
the universal trilinear scalar coupling Aq of comparable magnitude. The CMSSM parameter 
tan/3 ~ 35 and mi/2 -C mo, where mi/2 denotes the universal gaugino mass. Values of niQ 
and \Aq\ of comparable magnitude have been discussed from a different perspective string 
inspired models in Ref. 
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The search for NLSP stop, especially in the region of nearly degenerate stop and LSP 



neutralino masses, is challenging and has been implemented by both LEP and Tevatron 12 



14l |. In this case, the stop two-body decay channels into a top quark and neutralino, or a 
bottom quark and chargino, and three-body decay channels ii W~^bxi, ii — )■ bi^u and 
ii —J- bi^u are all kinematically forbidden. The loop-induced stop two-body decay into a 
charm quark and a neutralino is generally considered to overwhelm the four-body channel 
ti — )■ £^h'{qq')bxi and tends to be the dominant NLSP stop decay mode 15|]. The most 



stringent mass limit on a light stop with decay into charm quark and LSP neutralino comes 
from the CDF search for events containing two jets and missing transverse energy, namely 



M^^ > 180 GeV [IJ]. However, the Tevatron is not sensitive to stop searches if the stop and 
LSP neutralino mass difference is below 40 GeV. Thus the Tevatron bound does not cover 
the coannihilation region above the LEP limit of M^^ 100 GeV 
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Two alternative search methods have been investigated to detect a light stop. One of 
them takes advantage of the Majorana fermion feature of gluino and considers gluino pair 
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production followed by gluino decay into on-shell stop and top quark [l6|, [iTJ. The pair 
production of gluinos leads to events containing a pair of same-sign top quarks plus two same- 
sign stops. The benefit of this search is the anomalous same-sign dilepton signature arising 
from the same-sign top quarks leptonic decay, with negligible SM backgrounds. However, 
this method strongly depends on the production of relatively light gluino. From the well- 
known relation between gaugino masses at low energy in mSUGRA/CMSSM, namely M3 : 
M2 : Ml ^ Q : 2 : 1, which follows from the assumption of universal gaugino masses at 
high scale, the gluino masses have to be above at least 600 — 700 GeV for the coannihilation 
region of NLSP stop and LSP neutralino with M^^ > 100 GeV. This leads to small production 
cross sections, and so it is safe to assume that this scenario could not produce a significant 
amount of gluino pair events with same-sign dileptons at this early stage of LHC and evades 



the current bound especially for large values of mo 18|. The other proposed method is to 
consider stop pair production associated with a hard jet 19|]. In the coannihilation region, 
there will be minimal hadronic activity associated with the stop decay and therefore this 
channel would effectively lead to events with a hard jet and larg e missing energy. Such a 



signature has been proposed to explore large extra dimensions 20|, search for relatively light 



gluinos at Tevatron 2l|], and nearly degenerate gaugino pair production 22 1 



Recently, ATLAS reported results on an inclusive search for new phenomena in an event 
sample with monojet and large missing transverse momentum in the final state correspond- 
ing to an integrated luminosity of 1 fb~^ Good agreement was observed between the 
number of events in the data and the SM predictions. The results can be translated into im- 
proved limits on the stop-neutralino coannihilation scenario in terms of stop pair production 
associated with a hard jet. Also, the ATLAS and CMS experiments at ^/s = 7 TeV LHC 
have presented their results for events containing jets and missing transverse momentum for 



low-energy supersymmetry in 2011 (24,l25|. In mSUGRA/CMSSM with tan/3 = 10, Aq = 
and /i > 0, squarks and gluinos of equal mass are excluded for masses below 1 TeV or 
so. The lower limits on fundamental parameters of mSUGRA/CMSSM are mi/2 ~ 500 
GeV and mo ~ 3 TeV. However, the current "low statistics" with integrated luminosity 
up to 1 fb~^ encourages searches involving cascades initiated by gluino and the first two 
generation squarks, and consequently the hmits do not significantly depend on tan/3 and 



Aq parameters in mSUGRA/CMSSM 26|]. An analysis for scenarios with rich production 
of the third generation squarks induced by large values of mo and Aq in mSUGRA/CMSSM 
is still lacking. 



In this paper we study the LHC constraints on NLSP stop scenario in 6— r Yukawa unified 
mSUGRA/CMSSM using the monojet and multiple jets search results from the LHC. The 
paper is organized as follows. In section II we summarize the NLSP stop scenario with 
b — T Yukawa unification in mSUGRA/CMSSM and LSP neutralino (essentially bino) dark 
matter (for more details, see Ref. [10|). We also discuss the NLSP stop production modes 
and outline the selection cuts employed by the ATLAS collaboration. The results of the 
NLSP stop model constrained by the LHC data are presented in section HI. Our conclusions 
are summarized in section IV. 
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II. NLSP Stop in b-T Yukawa Unified mSUGRA/CMSSM and ATLAS Selection 

Requirements 



As mentioned earlier, the stop-bino coannihilation scenario requires the stop to be the 
NLSP in the sparticle spectrum, and to be nearly degenerate in mass with the bino LSP. 
The mass difference between the two should be 



Mi, 



^ < 20%. 



(1) 



This scenario can be realized in mSUGRA/CMSSM |27H29| because, compared to the other 



squarks, M^^ is reduced by contributions of the large top quark Yukawa coupling to the 
relevant renormalization group equations, as well as by mixing between the SU{2) doublet 
and singlet top squarks. To implement i/b = Ut at Mqut, sizable threshold correction to the 



bottom Yukawa coupling 6yb is required 10 



finite 



fx 



47r2 V 3 Ml 



gl Mg ^ y1 At 



8 M| 



tan /3, 



(2) 



where is the strong gauge coupling, Mg is the gluino mass, and Mi ^ (M^^ + M^J/2, 
M2 ^ (Mj2 + /i)/2. Also, the hierarchy Mg < Mj^,M^^;Mt^ < /i, is assumed. For 
/i > 0, to get the correct (negative) threshold correction, the contribution from the chargino 
loop (the second term of Eq. ([2])) should not only cancel the contribution from the gluino 
loop (the first term of Eq. ([2])), but it also must provide the correct negative sign for 5yh. 
This, it turns out, is achieved only with large values of mo and |^o|; with tan/3 ~ 35. 



The software package ISAJET 7.80 30|] was employed in Ref. 10| to scan over the relevant 
four parameters with /i > 0, as well as renormalization group evolution of gauge and Yukawa 
couplings and all soft parameters, and finally the computation of the physical masses of 
all particles. A large number of relevant phenomenological constraints such as BR{h — )■ 
57) jilj, BR{Bu ru) [il!], LEP II bound on the lightest Higgs and all the sparticle mass 
bounds {32! were also implemented. Note that recently the CDF collaboration has reported 
an excess in the rare decay 5^ — i- fi'^fi~ using 7 fb~^ of integrated luminosity and the 
measured central value of BR{B^ — j- fi'^fi~) is at least five times larger than the expected 



SM value 
1.14 fb-i 



33 



However, the combination of CMS and LHCb searches using 0.34 fb~^ and 
integrated luminosity has not confirmed this excess. Indeed, they provide a more 
stringent upper limit on the branching ratio, namely BR{B^ — n~^fi^) < 1.08 x 10~® at 



95% confidence level 3^. In our analysis we apply the upper limit from LHC. 



The degree of Yukawa unification is quantified by the parameter R 35 



R 



max{yb, Vr) 
mm{yb,yr) 



(3) 



In Ref. 10|, R is required to be < 1.1, so that b — T Yukawa unification holds at 10% level or 
better. It was shown that, to get good b — r Yukawa unification in mSUGRA/CMSSM, the 
parameters lie in the range of 5 TeV ^ mo ^ 20 TeV, 35 ^ tan /3 ^ 40, with |y4o/mo| ~ 2.3 
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and mi/2 <^ fno- The NLSP stop is quasi-degenerate in mass with LSP neutrahno, and in 
our study it will be lighter than the top quark. The gluino is about 6 times heavier than 
the LSP neutralino, while the remaining sfermions all have much heavier masses, namely 
greater than 5 TeV for M^^ , M^^ ^ , M^^ ^ and 10 TeV for the first two family squarks. 



Following the theoretical estimates in Refs. 15|, l27|, l36|] and experimental search assump- 



tion [iJ], we assume that the NLSP stop decays with 100% branching fraction into a charm 
quark and neutralino LSP. Also, the total decay widths of the models we consider are of 
order 10"^*^ GeV, which guarantees that the stops promptly decay in the detector (the decay 
length is too short to be displaced). 

For the NLSP stop scenario, the pure stop pair production is the leading source of stops. 
However, the small mass difference between the NLSP stop and LSP neutralino induces 
very soft charm jets and low missing energy from NLSP stop decay, which very likely evades 
the current LHC search bounds for multiple energetic jets. At best, only a tiny range of 
light NLSP stops could pass the relevant selection cuts because of the relatively large cross 
sections. It is therefore important to include the hard QCD emission at the matrix element 
level in order to provide a hard jet and large missing energy, and thus explore more stringent 
limits on the NLSP stop. In this scenario the heavier gluino essentially decays into on-shell 
stop plus top quark. The energetic particles from the top decay could compensate for the loss 
of events arising from the low cross section of heavy gluino production and NLSP stop decay 
followed by a soft jet. Based on these considerations, we generate hard scattering processes 
of gluino pair production and stop pair production together with the same processes with 
one extra jet at the matrix element level using Madgraph/Madevent js^] 

PP 99, iiil, fiiK (4) 
with the gluino decaying with 50% branching ratio into ttt and tii each. We use Pythia to 



include decays and parton showering and hadronization 38| , and PGS-4 to simulate the im- 



portant detector effects with ATLAS-like parameters [39[. We take care to correctly match 
(without double-counting) between matrix element and showering generation of additional 
jets. In Madgraph/Madevent running we implement MLM matching with P^^-ordered show- 
ers and the shower-i^^- scheme with Qcut = 100 GeV as described in Ref. 40 . The cross 
sections are normalized to the next-to- leading order output of Prospino 2.1 41 1. 

The ATLAS and CMS collaborations have reported data in terms of events containing 
jets and large missing transverse momentum in y/s = 7 TeV proton-proton collisions, cor- 



responding to an integrated luminosity of 1 fb~^ 2J, |25|]. Taking ATLAS as an example, 
the selected events are required to have a leading jet pt of at least 130 GeV and other 
multiple jets with px greater than 40 GeV. The quantity rrieff, the scalar sum of and 
the transverse momenta of the highest pt jets, must be more than 1 TeV in most of the 
selection modes. Apparently, these selection requirements are too stringent for our case with 
nearly degenerate NLSP stop and LSP neutralino. In our case, most of the jets come from 
the NLSP stop decay and are forced to be kinematically extremely soft. Thus, most of the 
events with a given stop mass would be eliminated by the selection cuts. We therefore expect 
that the upper limit on the excluded stop mass for nearly degenerate NLSP stop and LSP 
neutralino scenario from the multi-jets search would be lower as compared with the bounds 
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on the gluino and the first two family squarks, although the contribution of the additional 
jet in Eq. (jlj) would improve the situation. 

More importantly for this scenario, the ATLAS experiment has looked for monojet plus 
missing energy events with the same 1 fb^^ integrated luminosity 23|. They searched for 
one extremely hard jet, large missing energy and nothing else. No excess above the SM 
background expectation was observed. With more strict selection cuts and more data, new 
lower bounds on non-SM cross sections are obtained that are roughly 5 times more stringent 
than those from the 2010 data. This analysis has been used to make constraints on large 
extra dimensions 23|] and model- independent interactions of dark matter 42]. In our case, 
this search would be sensitive to stop pair production associated with one hard jet, followed 
by stop decay into a soft jet and missing energy. This can be employed, as we show here, 
to find useful constraints on the NLSP stop scenario with nearly degenerate stop and LSP 
neutralino masses. 

In the search for monojet plus large missing transverse momentum, the signal events are 
selected according to 3 different cut requirements, named "LP", "HP" and "VHP" 23|] as 
shown in Table [H The LP (HP) selection requires a jet with pr > 120 GeV {pt > 250 
GeV), |r/^'^*| < 2 in the final state, and Xt > 120 GeV (Xr > 220 GeV). Events with a 
second leading jet pt above 30 GeV (60 GeV) in the region |?7| < 4.5 are rejected. For the 
HP selection, the pt of the third leading jet must be less than 30 GeV, and an additional 
requirement on the azimuthal separation A0(jet,p^***) > 0.5 between the missing transverse 
momentum and the direction of the second leading jet is required. This cut is used to 
select events with the first and second jets going in roughly the same direction to reduce 
background from j{W -^)tv as stated in Refs. 19|, |43|. The VHP selection is defined with 
the same requirements as in the HP region, but with thresholds on the leading jet pt and 
increased up to 350 GeV and 300 GeV, respectively. Also, events with charged leptons 
are also rejected. The 95% C.L. upper limits on effective cross section (cross section times 
acceptance) for non-SM processes for signal region LP, HP, VHP are also shown in the 
last row of Table HI Following Ref. 4^ we apply a x acceptance > o"exp as the exclusion 
requirement for each model, where a is the relevant total cross section and the acceptance 
is the ratio of signal events after and before selection cuts which reflects the effects of 
experimental efficiency. Note that the upper limit on the second/third leading jet may 
remove some of our signal events once the hard leading jet is required. It is because the 
subleading jets in the roughly opposite direction may become correspondingly harder to 
balance the large px of the leading jet, although the sum of neutralino momenta leads to a 
large amount of missing energy. Nevertheless, we still apply this cut as we follow the ATLAS 
limits for (Xcxp, and it does help to efficiently suppress the W/Z + jets and tt backgrounds. 



III. ATLAS Constraints on NLSP Stop and Neutralino Dark Matter 

To study the LHC constraints on this class of models, we generate more than half a 
million models by scanning the parameter space. From these, 3705 models pass the various 
experimental constraints listed in section II, and they have acceptable Yukawa unification 
[R < 1.1) and NLSP stop. Note that in our analysis we focus on 983 of them corresponding 
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LP 


HP 


VHP 


Leading jet pt (GeV) 


> 120 


> 250 


> 350 


Second jets pt (GeV) 


< 30 


< 60 


< 60 


Third jets px (GeV) 




< 30 


< 30 


A0(pr-,j2) 




> 0.5 


> 0.5 


(GeV) 


> 120 


> 220 


> 300 


ATLAS cTexp (pb) 


1.7 


0.11 


0.035 



TABLE I. Summary of selection cuts and 95% C.L. upper limits on the effective cross section for 
non-SM processes for signal region LP, HP and VHP containing hnal states with monojet and 
missing transverse momentum with 1 fb~^ luminosity, following the ATLAS data analyses 231 ]. 



to < Mt = 173.3 GeV, which also happens to be favored by electroweak baryogenesis. 

In Fig. [1] we show ax acceptance vs. M^^ for the models with Yukawa unification and 
NLSP stop, after applying the requirements in the ATLAS monojet regions LP, HP and 
VHP. The constrained values of NLSP stop mass increase as the required pt of leading jet 
gets higher in the three different regions as shown in Tabled namely M^^ ~ 110 GeV for 
LP region, and M^^ ~ 160 GeV for the HP and VHP regions. It is because the emitted jet 
recoils against the two associated stops in the transverse direction to the beams, and thus 
its Pt is somewhat correlated with the relevant stop mass. Combining the exclusions from 
the three regions, an NLSP stop mass below 140 GeV is essentially excluded. 

Fig.[2]shows the exclusion plot in the M^o— M^-^ plane, with Yukawa unification and NLSP 
stop models (green circle). The top line corresponds to the kinematic bound of ti — j- cx? 
channel which is open below this line. The region below the bottom most straight line 
corresponds to the stop decay channel ii — >■ bW~^x^. In the region between these two lines, 
a stop decay into a charm quark and LSP neutralino is the unique channel for our study, since 
we assume that the 4-body channel is always highly suppressed. The coannihilation bounds 
from Eq. ([T]) are also displayed in this plot. One can see that Tevatron bound does not 
cover the coannihilation region. However, the ATLAS monojet search does make additional 
inroads beyond the Tevatron, denoted by red triangles. For coannihilation region, it is more 
sensitive as the mass difference between the NLSP stop and LSP neutralino decreases. For 
the region with M^^ > 140 GeV the monojet search loses its capability when the mass 
difference is larger than 20 GeV because the charm jets from stop decay in this case become 
harder and cannot pass the pt selection requirement for the non-leading jets. 

Besides the monojet channels, we also apply the ATLAS multi-jets search require- 
and show the excluded models (but not by monojet search) with black box in 



24 



ments 

Fig. [2j The requirement of the additional jet and heavier gluino also provides events with 
hard jet(s) and large missing energy that pass the multiple energetic jets search cuts. These 
excluded points gather in the region with M^^ — M^o > 20 GeV because of the induced 
relatively large pt of jets from stop decay. Based on these features we can clearly identify 
the region excluded by the LHC in Fig. |2j 

In Fig.[3l the excluded models are displayed in the mi/2 — 'tio plane with > and varying 
Aq and tan (3. We display all models that survive the low energy experiments listed in section 
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FIG. 1. cr X acceptance vs. Mj^ with liorizontal line as tlie 95% C.L. upper limits on effective non- 
SM processes cross section for signal region LP (top left), HP (top right), VHP (bottom). Green 
regions correspond to models with Yukawa unification (R < 1.1) and NLSP stop. 



II (grey color), models with good Yukawa unification and NLSP stop (green circle), and 
excluded models by the combined monojet and multi-jets searches (red triangle). One can 
see that the most stringent lower limit on tuq is around 3 TeV for tan /3 = 10, = 0, /i > 
case, from the LHC searches corresponding to comparable gluino and the first two family 
squarks masses 2J]. Our models with good Yukawa unification and NLSP stop correspond 
to mo > 8 TeV and the above study on the production of NLSP stop approach much larger 
values of mo, namely 8 TeV < mo < 16 TeV. A significant region of the parameter space is 
excluded. 

It is both interesting and important to see the implications of LHC data on direct and 
indirect dark matter detection in this class of Yukawa unified mSUGRA/CMSSM with NLSP 
stop. In Fig. IHwe display this by plotting the spin-independent and spin-dependent WIMP- 
nucleon scattering cross section asi (left panel) and asn (right panel) vs. M^o. A significant 
region around M^o ^ 100 GeV is excluded by LHC data, although it is allowed by CDMS-II, 
XENONIOO, SuperK and IceCube experiments. This excluded region even lies about one 
(six) order of magnitude below the expected XENON IT/SuperCDMS (IceCube DeepCore) 
bound for spin-independent (spin-dependent) cross section. 
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M,~, (GeV) 



FIG. 2. M^o vs. for models with Yukawa unification and NLSP stop (green circle), those 

excluded by ATLAS monojet regions (red triangle) and other excluded ones by ATLAS multi-jets 
regions (black box) in the framework of 6 — r Yukawa unified mSUGRA/CMSSM. The kinematic 
limits and coannihilation bounds are also displayed. The blue region refers to the excluded region 



by Tevatron 14]. 



IV. Summary 

Inspired by the recent LHC search for events containing monojet/multi-jets and large 
missing transverse momentum, corresponding to an integrated luminosity of 1 fb~^, we 
have explored its ramifications for mSUGRA/CMSSM models which display b — t Yukawa 
unification at 10% level or better, contain NLSP stop, and possess LSP neutralino dark 
matter. In this Yukawa unification framework, the approximate mass degeneracy between 
NLSP stop and LSP neutralino require relatively large values of mo (~ 8 — 20 TeV). This 
coannihilation scenario with NLSP stop decaying into a soft jet evades the previous Tevatron 
bound. In terms of the emission of a hard QCD jet associated with stop pair production, 
followed by stop decay into a soft charm quark and LSP neutralino, we find that the monojet 
search at ATLAS is sensitive to the region with small mass difference between NLSP stop 
and LSP neutralino. The excluded limit can reach 160 GeV for the NLSP stop mass in 
the coannihilation region, while NLSP stop mass below 140 GeV is essentially excluded. 
The analysis for the production of stops based on the above searches excludes a significant 
parameter region in mSUGRA/CMSSM, namely in the region 8 TeV ^ mo ^ 16 TeV. The 
LHC implications for spin- dependent and spin-independent LSP neutralino-nucleon cross 
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FIG. 3. mii2 vs. mo for models satisfying all low energy experiments (grey box), those with 
Yukawa unification and NLSP stop (green circle) and excluded ones by combined ATLAS monojet 
and multi-jets searches (red triangle) in the framework oih — r Yukawa unified mSUGRA/CMSSM. 
The most stringent bound on this plane from ATLAS is also displayed 24 1. 




FIG. 4. a SI (left panel) and ctsd (right panel) vs. M^o in the framework of 6 — r Yukawa 
unified mSUGRA/CMSSM. The excluded region is denoted in red. Current limits from CDMS-II, 
XENONIOO, SuperK and IceCube, and future projected sensitivities from XENONIT, SuperCDMS 
and IceCube DeepCore are also shown. 



sections are also explored. Regions of the parameter space, some lying well below the much 
anticipated future bounds from IceCube DeepCore, Xenon IT and SuperCDMS, are already 
excluded by utilizing the LHC data. 

Note added: After this paper was essentially finished, several papers related to the NLSP 
stop scenario have appeared |45l-l50 . 
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